Transcriptomic Analysis of Fiber Strength in Upland Cotton Chromosome Introgression Lines Carrying Different Gossypium barbadense Chromosomal Segments

Fiber strength is the key trait that determines fiber quality in cotton, and it is closely related to secondary cell wall synthesis. To understand the mechanism underlying fiber strength, we compared fiber transcriptomes from different G. barbadense chromosome introgression lines (CSILs) that had higher fiber strengths than their recipient, G. hirsutum acc. TM-1. A total of 18,288 differentially expressed genes (DEGs) were detected between CSIL-35431 and CSIL-31010, two CSILs with stronger fiber and TM-1 during secondary cell wall synthesis. Functional classification and enrichment analysis revealed that these DEGs were enriched for secondary cell wall biogenesis, glucuronoxylan biosynthesis, cellulose biosynthesis, sugar-mediated signaling pathways, and fatty acid biosynthesis. Pathway analysis showed that these DEGs participated in starch and sucrose metabolism (328 genes), glycolysis/gluconeogenesis (122 genes), phenylpropanoid biosynthesis (101 genes), and oxidative phosphorylation (87 genes), etc. Moreover, the expression of MYB- and NAC-type transcription factor genes were also dramatically different between the CSILs and TM-1. Being different to those of CSIL-31134, CSIL-35431 and CSIL-31010, there were many genes for fatty acid degradation and biosynthesis, and also for carbohydrate metabolism that were down-regulated in CSIL-35368. Metabolic pathway analysis in the CSILs showed that different pathways were changed, and some changes at the same developmental stage in some pathways. Our results extended our understanding that carbonhydrate metabolic pathway and secondary cell wall biosynthesis can affect the fiber strength and suggested more genes and/or pathways be related to complex fiber strength formation process.